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Abstract 



We describe a simple model for calculating the zero-temperature super- 
fluid density of Zn-doped YBa2Cu307_5 as a function of the fraction x of 
in-plane Cu atoms which are replaced by Zn. The basis of the calculation is 
a "Swiss cheese" picture of a single Cu02 layer, in which a substitutional Zn 
impurity creates a normal region of area Tr^ab around it as originally suggested 
by Nachumi et al. Here ^ab is the zero-temperature in-plane coherence length 
at X = 0. We use this picture to calculate the variation of the in-plane super- 
fluid density with x at temperature T = 0, using both a numerical approach 
and an analytical approximation. For 5 = 0.37, if we use the value ^ab = 18.3 
A, we find that the in-plane superfluid decreases with increasing x and van- 
ishes near Xc = 0.01 in the analytical approximation, and near Xc = 0.014 
in the numerical approach. Xc is quite sensitive to S^ab-, whose value is not 
widely agreed upon. The model also predicts a peak in the real part of the 
conductivity, Reo"e(<jj, x), at concentrations x ~ Xc, and low frequencies, and 
a variation of critical current density with x of the form Jdx) oc ns^eixY^^ 
near percolation, where ns,e{x) is the in-plane superfluid density. 

I. INTRODUCTION 

In typical s-wave superconductors, nonmagnetic impurities are not pair-breaking, and 
therefore have little effect on either the transition temperature or the superfluid density [|I| . 
By contrast, the substitution of a nonmagnetic impurity into a d-wave superconductor is 
expected to be pair-breaking, and, therefore, should have a significant effect on both. Since 
the Cu02-based high-temperature superconductors are believed to possess an order parame- 
ter with dx2-y2 symmetry 0], such an effect should be observable in these materials. Indeed, 
the substitution of a nonmagnetic impurity such as Zn (or even an insulating impurity such 
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as He, introduced by ion implantation) in a high-Tc material such as YBa2Cu307_5 is now 
well known to dramatically suppress the zero-temperature superfluid density and also 
to have a significant, though lesser, effect on the critical temperature. 

This striking behavior has provoked numerous theoretical studies of Zn impurities and 
other nonmagnetic disorder in the cuprate superconductors. For example, Annett et al 
and Hirschfeld and Goldenfeld |jTO| , have shown that the impurity scattering can convert the 
T-linear behavior of the low-temperature superfluid density to a behavior, in agreement 
with experiment. Other workers |TIHl3| have reflned the theory so that it now agrees quite 
well with experiment. Choi has used a t-matrix formalism to calculate the Zn impurity 
dependence of the superfluid density in YBa2Cu307_5 for various temperatures. Salluzzo 
et al [|l^] used a model of dirty d-wave superconductivity to flt the impurity-dependence 
of the superfluid density in a NdBa-based cuprate superconductor. Hettler and Hirschfeld 
T6| have shown how inhomogeneities in a d-wave order parameter can lead to both strong 



suppression of superfluid density and enhanced microwave conductivity. 

Several authors have treated the impurity problem by solving the Bogoliubov-de Gennes 
equations [0 for the pairing amplitude, which is position-dependent in a dirty d-wave or 
s-wave superconductor with a short coherence length. Ghosal et al used this approach 
to show that disorder in a 2D s-wave superconductor reduces the effective superfluid density, 
makes the superfluid inhomogeneous, and even breaks up the 2D layer into superconducting 
islands. More recently |18|, the same authors generalized this approach to d-wave super- 
conductors. Zhitomirsky and Walker |T9| showed how spatial fluctuations in the energy gap 
can reduce Tc in a d-wave superconductor. Haas et al ||2^ solved the Bogoliubov-de Gennes 
equations with a momentum- dependent pairing interaction to obtain an extended gapless 
region in the overdoped cuprates, which they connected to a rapid doping-induced decrease 
of both Tc and superfluid density. Franz et al , using the Bogoliubov-de Gennes approach 
for a spatially varying energy gap, found that the superfluid density goes rapidly to zero 
with increasing Zn concentration, while Tc decreases much more slowly with concentration. 



in reasonable agreement with experiment. Neto |^2[ has treated the Zn-doped cuprates as 
a collection of fluctuating stripes modeled as an array of coupled overdamped Josephson 
junctions pinned by impurities, and has used this model to estimate the doping-dependence 
of the transition temperature. Finally, Uemura has discussed the low-temperature su- 
perfluid density in Zn-doped cuprates from the point of view of the "universal correlation" 
between Tc and zero-temperature superfluid density in the underdoped cuprates. As noted 
by Emery and Kivelson such a connection would exist if were controlled by phase 
fluctuations rather than by the vanishing of the gap amplitude. 

Most of these models require a rather elaborate calculation to flnd the influence of non- 
magnetic impurities on the superfluid density. By contrast, the goal of this paper is to 
analyze the behavior to be expected of a very simple model of the Zn-doped cuprate super- 
conductors which can be worked out by intuitively appealing calculations. Although this 
model is difficult to derive from first principles, it may still be of value, because it can be 
used to make simple predictions for specific experiments, and to understand trends in those 
experiments. 



The model we study was first described by Nachumi et al and is sometimes known 
as the "Swiss cheese" model. In this model, the Zn atoms are assumed to replace the Cu 
atoms substitutionally within the Cu02 layer, which is believed to be the locus of super- 
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conductivity in the cuprate superconductors, and to disrupt the dx2_y2 order parameter, 
producing a region approximately of radius ^ab within which the layer is normal. (Here ^ab is 
the zero-temperature in-plane coherence length of YBa2Cu307_5.) This disruption is plau- 
sible, because ^ab is very small in the cuprate superconductors (probably on the order of 10 - 
20 A at low temperatures, as discussed further below). Indeed, some dramatic experimental 
evidence supporting both the Swiss cheese model and the d-wave order parameter picture 
has recently been presented by Pan et al [^. These workers showed, using scanning tun- 
neling microscopy, that the gap was significantly suppressed near the Zn atoms, and that 
the density of states of single-particle excitations near the Zn atoms exhibited a fourfold 
symmetry. 

In the present paper, we extend the approach of Nachumi et al by analyzing the Swiss 
cheese model within a simple percolation picture. Our approach differs from that of Nachumi 
et al by accounting for the connectivity of the superconducting part of the Cu02 layer. By 
contrast, in their approach, the superfluid density is directly proportional to the supercon- 
ducting area. The consequences of our approach are easily worked out if we use percolation 
theory and interpret the Cu02 layer as an inhomogeneous superconductor with spatially dis- 
tinct superconducting and normal region. The resulting variation of superfluid density with 
Zn concentration seems to agree somewhat better with experiment than does the simplest 
version of the Swiss-cheese model. 

The remainder of this paper is organized as follows. In the next section, we describe our 
percolation model in detail. The numerical results following from this model are presented 
in Section III. Finally, in Section IV, we discuss our numerical results, briefly compare them 
to experiment, and analyze the limitations and predictions of the model. 



II. MODEL 
A. Geometrical Assumptions 

Fig. 1 shows our picture of a single Cu02 layer in Zn-doped YBa2Cu307_5, in which 
a fraction x of the Cu atoms are replaced at random by Zn. Within the "Swiss cheese" 
model |2^, each Zn atom completely suppresses superconductivity in a region of area vr^^^ 
around it. Thus, as x increases, a smaller and smaller areal fraction of the Cu02 plane is 
superconducting. Eventually, for a large enough x, the remaining superconducting region is 
disconnected and the effective superfluid density must vanish. 

There are two parameters in our geometrical model: the fraction x of in-plane Cu atoms 
which are replaced by Zn, and the quantity ^ab/ci, which is the ratio of the zero-temperature 
in-plane coherence length ^ab (evaluated at x = 0) to the Cu-Cu distance a. Neither pa- 
rameter is trivial to determine experimentally. For example, x is not necessarily equal to 
c, the overall fraction of Cu atoms which are replaced by Zn. In fact, it is believed that 



X ~ 3c/2, i. e., that all the Zn dopant atoms go into the Cu02 planes [^]. Obviously, 
the connection between x and c is important in making quantitative comparisons between 
theory and experiment. Also, the magnitude of ^ab is not universally agreed upon for any 



value of S. Nachumi et al |^, for example, quote ^ab = 18. 3A at c = and S = 0.37, 



extrapolating from c = 0.01. Since a = 3.86A for this value of 6 |]2^, this gives ^ab/a = 4.74. 
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We shall later discuss the effects of assumptions about ^ab on the predictions obtained from 
this model. 



B. Assumptions about as and 



To estimate the superfiuid density as a function of x, we assume that the region more 
than a distance C,ab from any Zn atom has a complex conductivity (Js{uj) characteristic of a 
superconductor. We obtain this conductivity by using the first London equation, as applied 
to a Cu02 plane at x = 0: 

where is the supercurrent density and E is the electric field. In gaussian units, 
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where Xab is the in-plane penetration depth (evaluated at a; = 0), and m* = 2me is the mass 
and q = 2|e| the charge of a Cooper pair. Fourier transforming gives 3s{uj) = as{uj)'E{uj), 
where 

as{uj) = — (3) 

00 

with 

m* 4:TTXlh 

We assume that the region within a distance ^ab of a Zn atom is simply a normal metal with 
frequency-independent conductivity o"7v- 

Clearly, this simple model can only be appropriate for u <^ ljq, where o^o = A/aj^f. 
At higher frequencies, there will be effects arising from the superconducting energy gap, 
which will produce additional structure in as{uj). Also, this inductive form also implies 
no absorption at low frequencies in the limit x ^ 0. In reality, many of the cuprate 
superconductors do exhibit low-frequency absorption, even at low temperatures [^. This 
extra absorption can readily be included within the present model by adding to as a parallel 
normal conductivity. We have carried out a few calculations including this extra term; as 
discussed below, they do not affect the x-dependent superfiuid density. 



C. Definition of Eff"ective Superfiuid Density ns,e{x) 

We may infer the effective superfiuid density ns^eix) from an effective London equation 

A.flHE) (5) 
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which relates the space-averaged current density (J) to the space-averaged electric field (E). 
Ae(x) is related to ns^e{x) through a concentration-dependent generalization of eq. (2): 



47r 2 m* 



In the frequency domain, eqs. ^ and |^ may be written 



(JH) = ^^(EH) = aeM(E(..)), (7) 
where ae{uj,x) is the effective complex conductivity. Therefore, 

TTl* 

ns,e{^) = -^Lim^_,oUjImae{uj,x), (8) 

and, in order to compute nsA^)y need to calculate cre{uj, x) at sufficiently low frequencies. 

The ratio ns,e{^)/'^s,e{^) ^^^^ t)e calculated more directly by using a homogeneity 
property of (Xe ||3^. In the present context, this property states that 

Cre{<JS, <yN, X) = fiae{crs/fi, CTN/fJ', x) (9) 

Here ae{<Js,o'N,x) is the effective complex conductivity of a two-component system made 
up of constituents with conductivities as and and in-plane Zn concentration x in a 
particular geometry, and fi is an arbitrary constant. In effect, homogeneity means that 
(for fixed geometry and fixed x) if the conductivities of each constituent are multiplied by 
a certain factor the effective conductivity is multiplied by that same factor In 
particular, if // = cr^, we get 

0-e{<^S,(^N,x) = (JsO-e{l,aN/o-s,x). (lO) 

Hence, for any x, 

ns,e{^) Lim^^o^^Imcre{crs{uj),aN,x) 



ns,e(0) Limu;^oUjImae{crs{uj),aN,0) 



(11) 



ae(l,0,0)- ^ ^ 

Here we have used the fact that nsA^)y ^^e superfiuid density at x = 0, is given by ns^e = 
m*A/q'^, and also the fact that Lim^^oC"iv/cs(i-^) = 0. 



D. Two Methods for Calculating ns,e{^) 

We have evaluated ns^x) [or equivalently, ae{u!,x)] by two complementary methods. 
The first method is a finite-element technique, and should become very accurate for a suffi- 
ciently large numerical sample, given the geometrical assumptions of the model. We consider 
an L X L region of a Cu02 plane, assuming that the Cu ions sit on a square lattice of lattice 
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constant a, and choosing L = N^a where A^^,' is an integer. Then, using a random num- 
ber generator, we replace a fraction x of the Cu ions by Zn. We then assign conductances 
(Js{uj) or (Tat to each bond connecting nearest neighbor Cu sites, by the following rule. If 
the entire length of the bond is at least a distance ^ab from any Zn site, it is assigned a con- 
ductance crs{uj); otherwise, it is assigned (Jat. Next, we calculate ae{uj,x) using a standard 
technique known as the Y-A transformation |3^, which is very efficient and accurate for 
two-dimensional impedance networks. This technique is approximate only insofar as it is 
limited to finite networks. We call this method the random conductance approach (RCA). 

Our second method for calculating (Te{uj,x) is a simple analytical approximation, the 
effective medium approximation (EMA) pT|-p^. For the assumed geometry, this approxi- 
mation is implemented as follows. First, we calculate p, the areal fraction occupied by the 
superconducting material, as a function of x, using the relation 

p = Ltms-.ooll-^j , (13) 

where Nzn is the number of Zn ions in a plane of area S. This limit is readily evaluated 
using the relation x = NznCi'^/S, with the result 



p = exp -x^ . (14) 



a2 



Given p, one can calculate cJe using the standard EMA relation for this two-dimensional 
system, namely, 

p^^^ + {1 - p)^^^ = 0. (15) 

The physically correct solution to this quadratic equation is the one which satisfies ImcXe > 
0. Of course, ideally, the EMA is suited for application to a geometry in which the N 
and S components are distributed in a symmetrical fashion, whereas here, because of the 
assumption of excluded area, the N and S are not distributed symmetrically. Nonetheless, 
it is a reasonable approximation even for this nonsymmetric geometry. 

Once ae{uj) has been calculated at a given x by either of the methods described above, 
the effective superfluid density can be inferred from the relation (^. We can also compute 
the ratio ns,e{x)/ns,e{0) directly, using eq. ([T^). That is, for a given x, we decide which are 
the superconducting bonds, using the method described earlier in this section, and assign 
these bonds conductances of unity. The remaining normal bonds are assigned conductances 
of zero. The effective conductance is then evaluated using either the EMA or the RCA, and 
the ratio ^^s,e(a^)/''^5,e(0) is then obtained using eq. (|I^). 



III. RESULTS 

In carrying out the RCA calculations, we have generally used a square networks of size 
ranging from 160 x 160 to 500 x 500 Cu sites. In order to reduce statistical fluctuations, 
the RCA results quoted below are averaged over a number of different realizations of the 
disorder for each uj and x. These realizations are "correlated" in the sense that, for each 
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realization, we increase x by adding a few impurities to the configuration of the previous 
X value. In evaluating the EMA conductivities, no such averaging is necessary, since the 
approximation can be evaluated analytically in the limit of a very large sample. 

Our results for ^5,e(^)/^5,e(0) are shown in Fig. 2 for several choices of the parameter 
^afe/a. The full, dashed, dot-dashed, and dotted curves are the EMA results. The diamonds 
are the results of the RCA calculations at ^ab/d = 4.74, using eq. (12), for two correlated 
realizations of a 500 x 500 lattice. There are several striking features. First of all, in 
both methods ns^e{x)/ns,e{^) monotonically decreases with increasing x, reaching zero at a 
percolation threshold Xc which depends on ^ab/d- Moreover, Xc increases with decreasing 
C,ab/ CL. This relationship is not surprising: a larger ^ab/d implies that a larger area is converted 
from superconducting to normal by a single Zn impurity, within this model. Also, the EMA 
and RCA agree quite well over most of the concentration range, except near the percolation 
threshold. The EMA percolation threshold always corresponds to p = 0.5 |Q. The RCA 
threshold is higher, because the superconducting fraction remains connected to a somewhat 
higher fraction of Zn than is predicted by the EMA. 

In Figs. 3 and 4, we show Imcre(ti;, x)/crAr and Re(Te(ci;, x)/(j7v as a function of x at = 
O.OOltuo, where c^o = A/a^ and A is defined in eq. (H). We use using parameters thought 
to be appropriate for the case 5 = 0.37, namely, Cab/ a = 4.74 |2^. In both Figures, the 
dashed line represents the EMA; the diamonds denote the RCA (average over 20 correlated 
realizations of a 160 x 160 lattice) and the the open squares are the RCA, averaged over 
2 realizations of a 500 x 500 lattice. Fig. 3 shows that lmae{uJ,x)/aN closely mirrors the 
behavior of ns,e{x)/ns,e{0) shown in Fig. 2, as expected from the homogeneity relations (11) 
and (12). The apparent percolation threshold is again somewhat larger {xc ~ 0.014) in the 
RCA than in the EMA {xc ~ 0.01). For values of x well below Xc, lmae{uj = O.OOli-jQ, x)/aN, 
like ns^e{x) /ns^e{0), decreases linearly with x. This linear behavior is well known in the 
EMA 0. 

Fig. 4 shows that Iieae{u),x) has a strong peak occurring near the percolation threshold 
[xc ~ 0.01 in the EMA and ~ 0.014 within the RCA for our choice of ^ab/o)- The physical 
origin of this peak is discussed in the next section. We believe that the greater half-width 
of the RCA results arises because the calculation is carried out for a finite sample, and also 
because the average is taken over a number of different realizations, each of which has a 
slightly different percolation threshold for these finite samples. In support of this picture, 
note that the RCA half-width is smaller for the larger (500 x 500) samples, for which these 
fluctuations should be smaller. The slight random fluctuations in these curves as a function 
of X are, we beheve, also due to these finite-size effects. By contrast, the EMA results 
are obtained analytically for an effectively infinite sample; so these finite-size fiuctuations 
play no role. Note that, even for a 500 x 500 sample, there are still only 375 impurities at 
X = 0.015; so these fiuctuations are still significant in the RCA. 

In Fig. 5, we show Reae{uj,x)/aiy for u = O.OOOltUo, once again using ^ab/d = 4.74, as 
calculated in the EMA and in the RCA for two different sample sizes. The EMA peak 
is clearly stronger for uj = O.OOOIcuq than uj = O.OOIcjq- The RCA peak is also stronger, 
especially for the larger sample where the peak height is not washed out by finite size 
effects. The apparent double peak for the 500 x 500 sample is, we believe, just the result of 
finite-size fiuctuations, which are magnified at such low frequencies where the conductivity 
contrast \crs/cr]\f \ is larger. 
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Finally, we have calculated ae{uj^x)/(jN within the EMA for the same model as before, 
again with ^ab/O' = 4.74, but a slightly different choice for iJs{oj). Namely, we write 

o-5(t^) = — + a^, (16) 

UJ 

where we arbitrarily choose cr^ = ctn- The motivation behind this choice is that, in a 
material such as YBa2Cu307_5, which is believed to have an order parameter with dx'^-y^ 
symmetry, the in-plane conductivity is expected to have a quasiparticle contribution, even at 
very low temperatures. The additional term cr^ is a crude way of modeling this contribution. 
The resulting Heaeiuj.x) is plotted versus x aX oj = O.OOltuo, both for the model (0) and 
for a purely inductive as = iA/uj. Clearly, the peak near Xc is little affected by the change, 
but there is now, in addition, some extra contribution to Recre(ti;,x) for x < Xc- Although 
this extra contribution to as does slightly change Recre(ci;, x) at finite uj, the homogeneity 
relations (11) and (12) imply that it has no effect on ns,e{x), either in the EMA or in the 
RCA. 



IV. DISCUSSION 

It is of interest to compare our calculated values of Xc with the experimentally observed 
critical concentration for YBa2(Cui_cZnc)307_5 at 5 = 0.37. The critical value of c is 
between 0.02 and 0.03 ||7|,|3^-p^ (probably closer to 0.03, judging from in-plane resistivity 
measurements |^), as compared to our calculated critical value Xc ~ 0.014 for this 6 [cf. 
Fig. 2]. But a direct comparison of these values is difficult: the in-plane Zn concentration 
may differ from c. If in fact x ~ 3c/2, as mentioned above, then our model seems to 
underestimate Xc, the underestimate is less in the RCA than in the more approximate EMA. 
But this apparent discrepancy in Xc may not be very significant, for two reasons. First, Xc 
is quite sensitive to the parameter ^ab/d'- a change of ^ab only from 4.74 to 3 increases Xc 
by more than a factor of 2 in the EMA, and presumably by a similar amount in the RCA. 
As already mentioned, this parameter is not known experimentally with great precision. 
Secondly, our model assumes that the Cu02 layer undergoes an discontinuous change from 
S to N character at a distance ^ab from a Zn impurity. It would be more realistic to assume 
a gradual change. Although such a transition region might be difficult to include in the 
present model, its qualitative effect would probably be to reduce the effective radius of the 
normal region to a value smaller than ^ab- Such a reduction would, according to our results. 



increase Xc, making it closer to the reported value |37 



It is of interest to compare the present model to several others in the literature. In 



the original model of Nachumi et al |25], the superfluid density is proportional to the su- 
perconducting areal fraction, i. e., to the fraction of the Cu02 plane which is more than a 
distance ^ab from any Zn impurity. The present model improves on the implementation of 
this Swiss cheese model by taking into account effects of connectivity of the superconduct- 
ing parts. This treatment is consistent with the generalized London equation [eq. (5)]. In 



another class of models (see, e. g., Refs. IT3J2J) ns,e{x) is calculated using the Bogoliubov- 
de Gennes equation with some form of static disorder. This model, like the Swiss-cheese 
model, also leads to a spatially varying gap function or pairing potential. Although this 
model seems quite different from ours, we speculate that both actually have much of the 
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same physical content. The difference is that in the models of Refs. [0,^, the disorder is 
on an atomic scale and treated microscopically, whereas in the present case the disorder is 
on a greater length scale and treated within continuum electrodynamics. Nonetheless, both 
lead to surprisingly similar behavior (superfluid density diminishing rapidly with increasing 
X and vanishing at a critical concentration). We do not know how to connect the two on a 
more rigorous basis, however. 

A nontrivial prediction of the present model is the existence of a peak in Iieae{uj) near Xc- 
If such a peak were observed in experiments, it could be viewed as strong evidence in favor of 
the present model. This "percolation peak" is reminiscent of the well-known fluctuation peak 
in Recr(u;) near Tc in a homogeneous superconductor, but there are signiflcant differences in 
the underlying physics. The fluctuation peak results from slow motions of superconducting 
fluctuations with size of order ^ab(T). The spatial extent of these fluctuations diverges near 
Tc, as does the time scale on which they move; this diverging time scale leads to a strong 
low-frequency peak in the conductivity. By contrast the percolation peak involves motion 
of charge carriers through a static disordered structure; the frequency-dependence arises 
from the slow motion of these charges through the weakly connected structure (but static 
structure of superconductor that exists, in our model, near the pc (or Xc)- The characteristic 
length for this structure is the so-called percolation correlation length C,p, which diverges 
near pc on both sides of the percolation transition |Q . 

The existence of this percolation peak can be understood from the following crude argu- 
ment. On the normal metal side of the percolation threshold, the static d. c. conductivity 
is flnite, but should diverge as p pc- This divergence is due to the ever-larger super- 
conducting regions, which come closer and closer to shorting out the voltage across the 
normal metal as Pc is approached. This divergence should disappear as the frequency is 
increased, because the contrast in conductivities between superconducting and normal re- 
gions will become much smaller in absolute value (that is, the ratio Ictn/o'sI should approach 
unity). Therefore, at a flxed p just above pc (that is, on the normal metal side), Reae{uj,p), 
should show a strong maximum at uj = 0, and this maximum should diverge as p ^ pc- 
Since ae{uj,p) should vary continuously across the percolation threshold, except at exactly 
= 0, we expect the same behavior on the superconducting side of the threshold. This is 
qualitatively the behavior we see in both our EMA and RCA results. 

The behavior of this "critical peak" can be estimated from the standard scaling theory 
of the percolation threshold (see, for example, Ref. [Q). This scaling approach is based 
on the assumption that the percolation threshold is a critical point described by a single 
diverging length ^p] if that assumption is correct, then the scaling approach is believed to 
be exact. As applied to the present model, and using the homogeneity relations (11) and 
(12). this scaling approach dictates that cxe have the form 

ae(l, ar^/as, x) = as\^x\'F (^p) , (17) 

where Ax = x — Xc-, and s and t are standard percolation exponents. The function F{z) of 
the complex variable z has the limiting behaviors |^ 



F(2) ~ Ciz, 1^1 < 1; Ax > 0; (18) 
F(2)~C2, |^|<l;Ax<0; (19) 
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~C32*/(*+^),|2| >1, (20) 

where Ci, C2, and C3 are constants. This form imphes that precisely at x = Xc both the 
real and imaginary parts of de satisfy 

Reaeiuj, Xe), Imaeioo, x,) oc af+'^af'^'^ oc uj-'/'^'+'\ (21) 

In conventional 2D percolation theory, s = t ~ 1.30; so both Iieae{uj, Xc) and lmae{uj, Xc) 
would be proportional to c<j~^/^. The present "Swiss cheese" model may not, however, be in 
the same universality class as percolation on a lattice. Instead, it may fall into the class of 
certain continuum percolation models In the present case, the exponent t and s should 



have the values corresponding to another "Swiss cheese" model in 2D, in which the area 
between the holes is conducting while the holes are insulating. According to the discussion 
of Ref. [^, the Swiss cheese t and s are, in fact, unchanged from their lattice values of 1.3 in 



two dimensions. Therefore, one still expects that both Re(Te(c<J, Xc) and Im ae{uj,Xc) should 
have a cj"^/^ frequency dependence in this model. 

It is also of interest to consider the predictions of the present model regarding the critical 
current density Jc{x) in Zn-doped YBa2Cu307_5. In this case, we use the results of contin- 
uum percolation theory as applied to Jc{x) or equivalently, to Jcij))- two dimensions, 
for the present Swiss-cheese model, it is found that Jdj)) oc {p — PcY ■ The critical ex- 
ponent V is estimated as v = [u + l){d — 1) for Swiss cheese in d dimensions, where v is the 
percolation correlation length exponent defined by S,p oc \p — Pc\~'^ ■ Thus, in two dimensions 
V = u + 1. The conductivity exponent t is unchanged from its lattice value in c? = 2. Using 
the approximations t~z/ = 4/3in(i = 2 we obtain v ~ 7t/4. Equivalently, this may 
be written (for p near the percolation threshold) Jc{p) / Jc{^) ~ [ns,e{p)/'nsY^'^ ■ In terms 
of X, Jc{x)/Jc{0) oc [ns,eix)/ns,eiO)Y^'^- In summary, Jc{x) in Zn-doped YBa2Cu307_5 is 
predicted to go to zero near the percolation threshold even more rapidly than does ns,eix)- 
It would be of interest if this prediction could be tested experimentally. 

The present model does not give any information directly about the superconducting 
transition temperature Tc(x). Some experiments on Zn-doped YBa2Cu307_5 indicate that 
Tc(x) decreases much more slowly with x than does ns^e{x, T = 0) (see, for example, Ref. 0). 
Although we have no quantitative model for this behavior, a simple qualitative argument 
does suggest a explanation. Let us suppose that x < Xc, so that the superconducting portion 
of the material is connected throughout the plane. In that case, if thermal fluctuation effects 
can be neglected, then Tc(x) (as measured by a vanishing d. c. resistivity) should be indepen- 
dent of X for X < Xc, and zero for x > Xc- In this limit, clearly Tc(x)/Tc(0) > ns^e{.x) /ns^eif^)- 
In reality, thermal fluctuations should reduce Tc(x) below its expected value in the absence 
of fluctuations, but this inequality should still hold, as observed experimentally. 

In summary, we have presented a simple model for the variation of superfiuid density 
with in-plane Zn concentration x in Zn-doped YBa2Cu307_5. The model is based on the 



"Swiss cheese" picture suggested by Nachumi et al p5[. Using this model, we calculate 
the ratio ris^e{.x) /ris^eif^) for Zn-doped YBa2Cu307_5 as a function of the ratio ^ab/d, using 
two different approximations. The model predicts a critical Zn concentration above which 
ns^e{x) is zero, in rough agreement with experiment for Zn-doped YBa2Cu307_5 at 6 = 
0.37. The model also predicts a low-frequency peak in Reae{uj) near Xc, the existence of 
which has apparently not been checked experimentally. Finally, the model suggests that 
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the low-temperature critical current ratio Jc{x)/Jc{0) ~ [ns,e{x)/ns,eiO)Y^'^ near x = Xc- 
The observation of these features would give additional experimental support to the present 
simple model. 
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FIGURES 
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FIG. 1. Schematic diagram of model used to calculate effective superfluid density of Zn-doped 
YBa2Cu307_5. Shown is a sketch of a section of a Cu02 plane. The filled circles represent Cu 
ions. The open circles represent Zn ions, which have substitutionally replaced the Cu ions. O ions 
(not shown) are situated in the middle of each bond. The resistor-like and wavy lines represent 
superconducting and normal bonds in our model. The large circles, of radius ^ab, represent the 
regions which are assumed to have been driven normal by the Zn ions. The two diagonal lines 
represent the leads used in the application of the Y-A transformation. 
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FIG. 2. Relative zero-temperature superfiuid density ns,e{x)/ns,e{^)-, plotted as a function 
of in-plane Zn atomic concentration x for various values of the ratio ^ab/o,, where S^ah is the 
zero-temperature in-plane coherence length and a the lattice constant (Cu-Cu spacing). Dot- 
ted, solid, dot-dashed and dashed curves correspond to EMA calculations carried out within the 
"Swiss cheese" model as described in the text, for ^ab/o- = 3, 4.74, 6 and 10 respectively. Filled 
diamond symbols are results for ^ab/o- = 4.74 but calculated in the RCA. Line segments connecting 
the diamonds are guides to the eye. Note the difference in the percolation threshold: Xc is greater 
for the RCA than for the EMA. 
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FIG. 3. lm(Te{ijJ,x) / (Tn versus in-plane Zn atomic concentration x for u = O.OOlwo, where 
(Jo = A/aN-, assuming ^ab/o, = 4.74. We use as{oj) = iA/u, = const. Filled diamonds are RCA 
calculations for a 160 x 160 lattice, averaged over twenty correlated realizations; open squares are 
RCA for two correlated realizations of a 500 x 500 lattice). Solid line through these points is a 
guide to the eye. Dashed line: effective- medium approximation (EMA). 
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FIG. 4. Same as Fig. 3 but for Reae{u:,x)/(7j\f. Solid and dotted lines are guides to the eye. 
The RCA calculations for the larger size are carried out only between x = 0.0095 and x = 0.02. 
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FIG. 5. Same as Fig. 4 but for to = O.OOOltJo; Cab/a = 4.74. In this case, the full diamonds 
represent twenty correlated realizations of the RCA for a 140 x 140 lattice, and the open squares 
are two realizations for a lattice of size 500 x 500. This latter calculation is carried out only between 
X = 0.01 and x = 0.02. Light solid and dotted lines just connect the calculated points. Dashed 
line with no data points: EMA. 
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FIG. 6. Reae{uj,x)/a]\f as calculated in the EMA for two cases. Dotted curve: (7s{uj) = iA/uj, 
UN = const.; dashed curve: (Js{^) = iA/uj + cttv, cat = const. Units are same as in Figs. 2-5. 
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